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Abstract. The F-GAMMA program is a coordinated effort of several observing facilities to understand the 
AGN/blazar phenomenon via a multi-frequency monitoring approach, especially in the era of Fermi-GST. Some 
60 prominent sources are monitored monthly with the Effelsberg 100-m telescope, the IRAM 30-m telescope and 
more frequently but in a less uniform fashion at the APEX 12-m telescope, covering from 2.64 to 345 GHz. The pro- 
gram has been running since January 2007 and here some of its findings are summarized, (a) There are two major 
variability patterns that the spectra of sources follow, one spectral-evolution-dominated and one achromatic, (b) 
The FSRQs show higher brightness temperatures indicative of larger Doppler factors at play and (c) a statistically 
significant radio-7-ray correlation has been found with a method recently suggested by Pavlidou et al. (in prep.). 
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According to our current and rather complete understand- 
ing, in terms of system configuration, of the magnificent 
phenomenon of Active Galactic Nuclei (AGNs) , the incon- 
ceivable amounts of energy radiated by these systems are 
generated by a super-massive black hole living in the nu- 
clear region of the host galaxy. The even more extended 
variety of AGN phenomenologies are attributed to the 
same system seen at different viewing an gles by the admit- 
tedly elegant unified model, reviewed bv lUrrv fc Padovani 
(I1995I) . 31azars, as the most dramatic manifestation of the 
AGN phenomenon, display extreme observational charac- 
teristics due to their close to our line-of-sight orientation 
and the consequent boosting effects. Their double humped 
Spectral Energy Distribution (SED) is explained by means 
of two components. The low energy peaked one is ex- 
plained as a synchrotron component whereas the high en- 
ergy peaked one is assumed to be the result of inverse 
Compton processes either on external seed photons (from 
the disc, the torus or the broad line region) or synchrotron 
photons. What the exact mechanism is, is still unclear. 
Furthermore, what exactly causes the variability at vast 
timescale ranges is also unknown and so are numbers of 
fundamental questions. 

Among the several incredible opportunities offered by 
Fermi-GST is the densely sampled 7-ray light curves and 
spectra it provides and which for the first time allow re- 
ally simultaneous multi-energy studies. The F-GAMMA 
program is utilizing this approach for understanding some 
of these fundamental questions. In particular, roughly 60 
sources have been monitored monthly since January 2007 
with the Effelsberg 100-m telescope, the IRAM 30-m tele- 
scope and the APEX 12-m telescope covering vast fre- 
quency range. The observations at different telescopes are 
coordinated within 1 week. Here a few findings are dis- 
cussed after a short introduction to the program. 



2. The F-GAMMA program 



It has been discussed elsewhere (jFuhrmann et al.l . 12007 : 
Angelakis et~all . l2009l) that the F-GAMMA program is 



the coordinated effort for a multi-frequency study of AGN 
astrophysics in the light of Fermi-GST gathering data 
since January 2007. The pivotal facilities are the Effelsberg 
100-m telescope covering the range between 2.64 and 
43.00 GHz in 8 frequency steps, the IRAM 30-m telescope 
covering 86, 142, 220 GHz and the Apex 12-m telescope 
operating at 345 GHz. The monitoring is done monthly 
on a sample of 60 sources. The observations at these fa- 
cilities are coordinated within roughly a week. Along with 
this effort there take place several other closely collaborat- 
ing programs such as the OVRO monitoring program at 
15 GHz, the IRAM 30-m polarization m onitoring program 
(Thum et al. and lAgudo et al. ( 2010l)). and the Perugia 

all . l2002t 



AIT blazar monitoring (jTosti et al.l . l2002t ICiprini et al 



I2OO8I ). Here we present Effelsberg data alone. The light 
curves produced and animated spectra are constantly up- 
dated at www.mpifr.de/div/vlbi/fgamma. 

3. The sample 

The main motivation for this work has been the multi- 
band approach to AGN physics and especially the blazar 
phenomenon, in particular the investigation of correla- 
tions between the radio and 7-ray emission. For this rea- 
son the first sample consisted of 61 blazars (32 FSRQs, 
23 BLLacs, 3 radio galaxies and 3 unclassified blazars) 
selected mostly on their likelihood of being Fermi-GST 
detectable (earlier EGRET detections) . The release of the 
LAT Bright AGN Source hst (LB AS, lAbdo et all . \200^ 
showed that 29 of the 61 sources (~ 47% ) were detected 
by LAT in the first three months of operation. In a major 
revision of the source sample, LBAS as well a s the Fermi 
LAT First Source Catalog, were consulted (jAbdo et al. 
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2010[) in order to build up a new source sample with max- 
imum Fermi-LAT detectability and reasonably fast and 
intense variability. Currently the observed sample includes 
a total of 65 sources from which 60 are observed monthly. 
From these, 36 are clarified as FSRQs, 17 as BL Lacs and 
9 as unclassified blazars (classification by iMassaro et al.l . 
l2005l[200& ,20091 Additionally, one source is classified as 
a Seyfert and three are radio galaxies. 

From the above discussion it is very clear that the F- 
GAMMA sample suffers severely from biases and there- 
fore it is statistically incomplete. Nevertheless, extracted 
generalizations may be tested by careful comparison with 
other statistically complete samples. In any case, any gen- 
eralization must be stated with caution. 



4. Spectral Variability 

Blazar variability has been attributed to several factors 
and accordingly certain models have been developed to 
interpret and quantitatively describe it. One could dis- 
tinguish between two model categories: (a) models which 
attribute the variability to mechanisms that impose spec- 
tral evol ution on the observed ev ents such as shock-in-jet 

1985ll 



models (Marscher fc Gear 



plasma shells (|Guetta et al 



2004f) and (b) models which 



explain the variability in terms of geometrical effects such 
as systematic changes in the beam orientation (e.g . light- 
house effect, ICamenzind fc Krockenberger , 1992) which 
could for instance cause a change in the Doppler factor and 
further induce achromatic changes in the observed radio 
spectrum. The expected bi-modality in the phenomeno- 
logical behavior of the observed radio spectra is indeed 
seen in the f -GAMMA data. 

Figure [1] shows the EfFelsberg light curves and radio 
spectra for two representative cases that resemble what 
would be expected from the previous classification. In the 
case of 0235-1-164 the variability is dominated by spec- 
tral evolution (hereafter type 1), evidence of a three-stage 
evolut ionary path (Conipton, synchrotron and adiabatic 
losses. iMarscher fc GeaA ll985l ). The case of 0814-^425 on 
the other hand is representative of achromatic variability 
indicating some geometrical effect (hereafter type 2). 

Interestingly, from the point of view of the variability 
pattern followed by the observed spectra, all the sources 
fall in only five classes which comprise modifications of 
these two basic behaviors. 

The fact that there must be a fundamental differ- 
ence in the mechanism causing the variability becomes 
evident also from the evolutionary paths followed by the 
turnover frequency and fiux density at that point {Sm, 
Vin)- In figure [5] are shown the evolutionary tracks for 
three cases of type 1 and three cases of type 2 after the 
subtraction of mean quiescent spectra (with S (x i>°' and 
a = —0.5). The former c ase seems to be described well by 
Marscher fc Geaii (|l985l ) whereas type 2, needs a different 
interpretation. In a forthcoming publication (Angelakis et 
al. in prep.) we use this approach to estimate source pa- 
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Fig. 1. Two prototype cases representing the two charac- 
teristic behaviors. In the case of 0235 + 164 (1"' and 2"^^ 
panels from the top) the variability is dominated spec- 
tral evolution. The case of 0814 -I- 425 is representative of 
achromatic variability indicating some geometrical effect. 



rameters, and investigate the presence of possible quasi- 
periodicities in type 2 sources 

From average spectra a low and high frequency spec- 
tral index can be calculated in order to examine whether 
there is a significant differentiation between BL Lacs and 
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Fig. 3. The distribution of high-frequency spectral indices 
for BL Lacs and FSRQs with the former centering around 
flatter values than the latter. 
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Fig. 2. The evolution of the peaks of convex spectra in 
the Sm — Vm space for three members of each of the two 
characteristic classes type 1 and type 2. In every case a qui- 
escent spectrum of —0.5 is assumed {S oc v"). Members of 
the same class follow qualitatively the same evolutionary 
paths. 

FSRQs. Figure [3] shows the spectral indices between 9 mm 
observed with the Effelsberg telescope and 3 and 2 mm ob- 
tained with the IRAM 30-m telescope. What seems to be 
the case is a tendency of the BL Lacs to center around flat- 
ter high-frequency spectral indices as compared to FSRQs 
(and is indicated by statistical tests). This may for in- 
stance be due to the fact that the BL Lacs systematically 
peak at higher frequencies or, alternatively, their flares 
expand before reaching lower frequencies, maybe due to 
lower power. This requires further investigations. 

5. Time series analysis 

The level of variability present at all frequencies and in 
all sources becomes obvious as excess variance and has 
been quantifled via a formal t^st for which a signifl- 
cance reference level of 99.9% has been adopted. For the 
first 2.5 years of observations (that is for the first target 
sample), tests show that almost all the sources ap- 
pear significantly variable at all wavelengths (91 % of the 
sources are significantly variable with the proportion of 
variable sources dropping towards higher frequencies due 
to larger measurement uncertainties at these frequencies) . 



Using the mean rms ((cr)) averaged over all sources at 
each frequency as a measure of the variability amplitude, 
one can clearly see a monotonically rising trend as a func- 
tion of frequency. That is expected from sources which are 
mostly spectral-evolution-dominated. 

In order to investigate the presence of characteris- 
tic timescales in the acquired light curv es, the fi r st or- 
der Structure Function has been used (IRutmanl . [1978; 
Paltani et aLl . ll997t[si^onetti et allll985[) . From the time 
series analysis applied to the first 2.5 years of data, the 
variability timescales range between 80 and 500 days. 
From the estimated timescales and on the basis of some 
fundamental assumptions (causality, a single emitting 
component), one can estimate the brightness temperatures 

Tb = 4.5 • 10^° AS (^:^^^y, with Tb in K, AS" the fiux 
density variations in Jy, A in cm, Dl the luminosity dis- 
tance in Mpc, T the characteristic timescale in days and 
z the redshift. Figure |4] shows the maximum brightness 
temperature distribution (calculated from the fastest time 
scales reliably detected) for FSRQs and BL Lacs. There is 
a clear and statistically significant separation between the 
two classes, with the FSRQs showing brightness temper- 
atures systematically higher than the BLLacs. A possi- 
ble explanation for this could be that the former undergo 
stronger Doppler boosting. 

6. Radio versus 7-ray fluxes 

One of the remaining fundamental questions in AGN as- 
trophysics is where and how high energy emission is pro- 
duced. Searches for correlations between radio and 7-ray 
luminosities is used to study the connection between the 
mechanisms producing them and has been a rather de- 
batable field. The situation becomes even more perplex- 
ing due to artifacts that may be introduced by selection 
biases, redshift dependencies, lack of simultaneity and so 
forth. 

For a sample of 29 i^-GAMMA sources from the first 
sample detected by Fermi-GST, simultaneous 7-ray and 
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Fig. 4. The brightness temperature distribution for 
FSRQs and BL Lacs separately. 

radio data have been used to investigate the presence of 
such a correlation. Figure [5] shows that there indeed ex- 
ists a correlation. A new Monte-Carlo method has been 
used to assess its significance and it shows that it is in- 
deed intrinsic. The details of the newly applied method 
are describe by Pavlidou et al. (in prep.). 
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Fig. 5. Flux- flux correlation plots with three-month aver- 
aged radio data at 4.85 GHz. 



7. Conclusions 

After 3.5 years of observing, the f -GAMMA program has 
produced a large volume of data which allows very detailed 
blazar studies. Some of them have been discussed here: 

— The spectrum variability is not erratic but follows only 
two behaviors: (a) one dominated by spectral evolu- 



tion and (b) one achromatic variability pattern resem- 
bling rather a geometrically induced variability. All the 
sources fall in these categories and their modifications. 
Differences in the variability properties of sources from 
different classes must be sought. 

— The high-frequency spectral indices imply a separa- 
tion between FSRQs and BL Lacs, possibly due to the 
BL Lacs peaking at higher frequencies. 

— A clear distinction in the brightness temperature dis- 
tributions for BLLacs and FSRQs may be due to 
higher Dopplcr factors at play in FSRQs. 

— A newly suggested Monte-Carlo method assesses sta- 
tistically significant radio and gamma-ray flux correla- 
tions. 
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